Astron. Nachr. / AN 326, No. 9, 864-866 (2005) /DOI 10. 1002/asna.2005 10428 



The correlations and anticorrelations in QPO data 



M. A. Abramowicz 1 , D. Barret 2 , M. Bursa 3 , J. Horak 3 , W. Kluzniak 4 , P. 
Rebusco 5 , G. Torok 6 

1 Department of Physics, Goteborg University, S 412 96 Goteborg, Sweden 

2 Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, 9 Avenue du Colonel Roche, 31028 Tolouse Cedex 04, France 

3 Astronomical Institute, Academy of Sciences, Bocm II 1401, 141 31 Praha 4, Czechia 

4 Max-Planck-Institute fiir Astrophysik, D-85741 Garching, Germany 

5 Zielona Gora University, Lubuska 2, 65-265 Zielona Gora, Poland 

6 Institute of Physics, Silesian University in Opava, Bezrucovo nam. 13, 746 01 Opava, Czechia 



Received; accepted; published online 



Abstract. Double peak kHz QPO frequencies in neutron star sources varies in time by a factor of hundreds Hz while in 
microquasar sources the frequencies are fixed and located at the line U2 = 1.5^i in the frequency-frequency plot. The crucial 
question in the theory of twin HFQPOs is whether or not those observed in neutron-star systems are essentially different 
from those observed in black holes. In black hole systems the twin HFQPOs are known to be in a 3:2 ratio for each source. 
At first sight, this seems not to be the case for neutron stars. For each individual neutron star, the upper and lower kHz QPO 
frequencies, vi and v\, are linearly correlated, — A v\ + B, with the slope A < 1.5, i.e., the frequencies definitely are 
not in a 1.5 ratio. In this contribution we show that when considered jointly on a frequency-frequency plot, the data for the 
twin kHz QPO frequencies in several (as opposed to one) neutron stars uniquely pick out a certain preferred frequency ratio 
that is equal to 1.5 for the six sources examined so far. 



Key words: QPOs - Neutron stars 

©2005 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 

1. Introduction 

It has been recognized some time ago (Bursa 2002) that the 
pairs of high-frequency QPOs observed in all neutron star 
sources can be fitted by a single linear relation 
v 2 = 0.89i/i +375 Hz , (1) 
where v 2 and v\ refers to the upper and lower observed 
kHz QPO frequency. The linear correlation between the 
pair frequencies becomes even more apparent if we restrict 
ourselves to individual sources, though the coefficients of 
the relation Q slightly differ from source to source. A 
particular example is the case of ScoX-1, where the ob- 
served frequencies are well fitted with the linear dependence 
v 2 = (3 /4) v\ + 450 Hz with the accuracy better than one 

per cent. 

lAbramowicz et all J2003h pointed out that the distribution 
of frequency ratios is strongly peaked near the 1.5 value, but 
made no statement about the actual correlation of frequen- 
cies. A concern has been raised, wh ether their result is a sta- 
tistical illusion iBelloni etall20 05). 
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In this contribution we examine the linear frequency- 
frequency correlations for six neutron-star sources: 
4U 1820-30, 4U 1728-34, 4U 0614+09, 4U 1608-52, 
4U 1735-44 and 4U 1636-53, and show that within errors 
the correlations are consistent with the following statement: 
the plots of the six v 2 = Av\ + B relations intersect in 
one point [N\,N2\, with A^i«600Hz, and the ratio of 
the frequencies in the intersection point is N2/Ni~1.5 to 
high accuracy. To demonstrate this, we plot the five pairs 
of coefficients, A and B, and show that they are linearly 
antic orrelated. 

2. Anticorrelation between the slope and shift 

We examine the six individual neutron star sources by fitting 
each of them with a linear formula 

v 2 =Av x +B, (2) 

where the coefficients A and B are referred to as the slope and 
shift, respectively. The resulting values of the slope and shift 
and of corresponding errors for each source are summarized 
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Table 1. Best linear fits and their errors for the frequency- 
frequency correlation. 



l.G 



Source 


A 


AA 


B[Hz] 


AB [Hz] 


4U 1636 


0.58 


0.03 


622 


27 


4U 1608 


0.76 


0.03 


456 


16 


4U 1820 


1.02 


0.09 


255 


66 


4U 1735 


0.61 


0.05 


593 


39 


4U 0614 


0.79 


0.02 


420 


5 


4U 1728 


0.98 


0.02 


330 


5 



in Table [2 and plotted in Figure [2 showing the slope-shift 
plane. The dependence A — A(B) strongly suggests that the 
two quantities are anticorrelated. The linear fit for the anti- 
correlation gives 



A = (1.50 ± 0.03) - (0.0016 ± 0.0001) Hz" 1 B 



(3) 



As it was mentioned above, this result is consistent with the 
statement that the six linear plots of the V2-V1 relations inter- 
sect in one point [N2, N\]. The frequencies N2 and Ni can 
be determined from the coefficients of the anticorrelation Q, 
from which we obtain 

N 2 = (940 ± 80) Hz , TYi = (625 ± 40) Hz . (4) 
The ratio of frequencies N2/N1 at the intersection point 
equals to 1.5 with the accuracy of 15%. Hence, it is remark- 
able that neutron stars QPOs pick up the same 3:2 ratio com- 
mon for the black-hole sources. We can rephrase the equation 
(|5Jl to the form 
3 B 



A = — — 

2 625 Hz 



(5) 



The data points of Figure 1 have been obtained through 
a shift-and-add technique, as described in iBarret et alj j2005) 
applied to the whole archival RXTE data available to date. 
The technique relies firstly on estimating the QPO frequency 
down to the shortest timescales permitted by the data stati- 
tistics, and secondly on identifying the QPOs for which the 
shift-and-add can be applied (generally the lower QPOs). 
Having a set of frequencies, the frequencies are then binned 
(typically with a bin of 10 to 20 Hz), and the individual 
PDS shifted to the mean frequency within the bin. The aver- 
aged PDS so-obtained is then searched for QPOs, and when 
both QPOs are present, their frequencies are obtained with 
a Lorentzian fit. A linear fit is then applied to the lower and 
upper QPO frequencies so obtained. This method has been 
applied to 4 sources shown in Figure 1 (4U1820; 4U1608, 
4U1735, 4U1636), whereas for 0614 and 1728, the linear fit 
has been applied to the frequencies measured with a shift and 
add technique applied to a continuous segment of observation 
(of typical duration close to the orbital period of the RXTE 
spacecraft). Both methods should yield consistent results: the 
smaller error bars the points of 0614 and 1728 can be ex- 
plained by a larger number of frequencies involved in the lin- 
ear fit. A more careful comparison of the results of the two 
methods (with estimates of the potential biases, such as only 
the narrower QPOs are detected, hence with small errors on 
the frequency determination) is clearly required and will be 
part of a forthcoming paper (Barret et al. 2005b in prepara- 
tion). 
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Fig. 1. The anticorrelation between the shift and slope and the 
effect of 1/M-scaling. The figure shows the result obtained 
from the analysis of the data. Each point in the plot corre- 
sponds to the particular source. The shift A and the slope 
B correspond to the best linear fit of the v\—v% correlation. 
Clearly, A and B are anticorrelated among the sources with 
the best fit (solid line) given by equation (|3}- In the spirit 
of 1/M scaling, it corresponds to a 'typical' neutron-star 
mass between the examined system. Dashed lines represent 
the best fit recalled to 0.92M and 1.25M. The value of the 
fit crosses vertical axis very close 1.5, which is in agreement 
with the idea that the excitation mechanism of QPOs is the 
3:2 resonance. 

3. A possible explanation 

It has been suggested that QPO arises from a resonant inter- 
action between the radial and vertical oscillation modes in 
relativistic accretion flows. The strongest possible resonance 
occurs at the radius, where the radial and vertical epicyclic 
frequencies are in the 3:2 ratio. 

In general, the frequency and the amplitude of non-linear 
oscillations are not independent. In the lowest order of ap- 
proximation the observed frequencies v differ from the eigen- 
frequencies vq of oscillators by corrections proportional to 
the squared amplitudes. We consider system having two os- 
cillation modes, whose eigenfrequencies are v\q and i^.o- 
The frequencies of non-linear oscillations may be written in 
the form 



(6) 



The frequency corrections Az/i and Av 2 are proportional to 
the squared amplitudes 



Avx = n\a 1 + K2CL2 1 ^ v 2 = Aid} + \ 2 o-2 



(1) 
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where k%, k%, \\ and A2 are constants depending on non- 
linearities in the system. 

Let us suppose that for some reason the two amplitudes a\ 
and a,2 are correlated. Thus, one may consider the amplitudes 
as functions of a single parameter s, 



fti = ai(s) , a 2 = a 2 (s) . 



(8) 



Such relations can be considered as a natural consequence 
of an interplay between the resonance excitation mechanism 
and the dissipation of the energy in the system. 

It follows that the frequencies of non-linear oscillations 
v\ and V2 are correlated as well. Up to the linear order in s, 
we obtain from equation (J6J 

vi = fi,o + ^1,0*13 + 0{s 2 ), (9) 



V2 = V 2 fi + V 2 fiF 2 S + 0(s 2 ) . 



(10) 



where the coefficients F\ and F2 are given in terms of con- 
stants «i, K2, Ai and A 2 of frequency corrections and of 
the derivatives dai/ds and da 2 /ds at the point s = 0. 

Isolating the parameter s from equations (|9j and d 1 01 we 
get the linear correlation between the observed frequencies 

V2=Avi+B, (11) 

where the slope and shift are respectively given as 

A=^Q, (12) 

B = 1/2,0(1 ~ Q) (13) 

and we define Q as Q = F 2 jF\ . 

The observed frequencies v\ and v 2 of systems with dif- 
ferent amplitude prescription (jSJl are correlated in a different 
way. Any particular value of Q leads to particular values of 
the slope and the shift. However, if the eigenfrequencies of 
the systems are similar, the slope and shift are necessarily an- 
ticorrelated. Solving the equations ill 2b and Jl 3I > for the pa- 
rameter Q, we obtain 

A= ^o_J_ B (14) 

From the resonance condition it follows that the eigenfre- 
quency ratio ^2,0/^1,0 is approximately 3/2. Therefore we ar- 
rive at 

A=l-—B. (15) 

2 Vl,Q 



are close to 3:2 ratio supports the idea that there is a simi- 
lar mechanism at work in both classes of sources. Moreover, 
if we extend the black-hole 1/M scaling law up to the fre- 
quency of the intersection point, we obtain a mass of order of 
one solar mass (assuming zero angular momentum for neu- 
tron stars), which provides an additional hint. 

Assuming that the 1/Af-scaling can be adopted also to 
the neutron-star QPOs, the fact that the individual positions 
of sources in the A—B plane do not strictly follow the anticor- 
relation line can be attributed to small differences in neutron- 
star masses. By scaling the 614 Hz frequency of equation (J3j, 
we find that the A-B is steeper or softer for more massive or 
less massive sources, respectively. This is demonstrated by 
the shaded region in Figure^ Under this assumption the de- 
viation in the masses of examined neutron stars should not be 
greater than ~ 20%. 
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4. Discussion and Conclusions 



In black-hole sources the observed QPO frequencies are fixed 
and always have the ratio 3/2. It has been recognized that 
their actual frequencies scales in versely with mass M as- 
sumin g a similar value of the spin jMcClintock & Rem illard 
I2005L [van der Klis], [Torok] 1 ). In neutron-star sources the 
frequencies are not fixed, but their distribution seems to clus- 
ter around a single line for each individual source. By linear 
fitting the observed data, we have found out that these lines 
intersect around a single point \N\,N^[, which have coor- 
dinates given by equation @. The fact that the frequencies 
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